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We study the number of lattice points in integer dilates of the rational poly- 
tope 

V = <.{xi,...,Xn) e K>o : ^ x^ak < 1 ^ 
I 1 
where ai, . . . ,a„ are positive integers. This polytope is closely related to the 
linear Diophantine problem of Frobenius: given relatively prime positive inte- 
gers ai,...,a„, find the largest value of t (the Frobenius number) such that 
mifli -I- • ■ • -I- mna„ — t has no solution in positive integers mi, . . . ,m„. This 
is equivalent to the problem of finding the largest dilate tV such that the facet 
{X]fc=i ^kOk ~ t} contains no lattice point. We present two methods for com- 
puting the Ehrhart quasipolynomials L{V,t) := #{tV n Z") and L{'P°,t) := 
#(f'P° n Z"). Within the computations a Dedekind- like finite Fourier sum ap- 
pears. We obtain a reciprocity law for these sums, generalizing a theorem of 
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Gesscl. As a corollary of our formulas, wc rcdorive the reciprocity law for Za- 
gier's higher-dimensional Dedekind sums. Finally, we find bounds for the Fourier- 
Dedekind sums and use them to give new bounds for the Frobenius number. 

Key Words: rational polytopes, lattice points, the linear diophantine 
problem of Frobenius, Ehrhart quasipolynomial, Dedekind sums 
2000 Mathematics Subject Classification: 11D04, 05A15, 11H06 

1. INTRODUCTION 

Let ai, . . . , a„ be positive integers, Z" C R" be the n-dimensional inte- 
ger lattice, and 

V=l^{xu...,Xn)€W: a;fe >0,^afca;fc < 1,| , (1) 
a rational polytopc with vertices 

(o,...,o),(^,o,...o),(o,^,o,...,o),...,(o,...,o,^) . 

For a positive integer t e N, let L{'P,t) be the number of lattice points 
in the dilated polytopc tP = {tx : x E V}. Denote further the relative 
interior of V by V° and the number of lattice points in tV° by L{P°,t). 
Then L{V°,t) and L{V,t) are quasipolynomials in t of degree n [11], i.e. 
expressions 

Cnit) f' +--- + Ci{t) t + Co(t), 

where cq, . . . ,c„ are periodic functions in t. In fact, if the a/c's are pair- 
wise relatively prime then ci , . . . c„ are constants, so only cq will show this 
periodic dependency on t. 

Let A = {ai, . . . , a„} be a set of relatively prime positive integers, and 

PaW = # . . . ,m„) e N" : ^ mfettfe = i| . (2) 

The function (t) can be described as the number of restricted partitions 
of t with parts in A, where we require that each part is used at least once. 
(We reserve the name pa for the enumeration function of those partitions 
which do not have this restriction.) Geometrically, p'j^{t) enumerates the 
lattice points on the skewed facet of V. Define /(ai,...,a„) to be the 
largest value of t for which 

p'At) = . 

In the 19th century, Frobenius inaugurated the study of /(ai, . . . , a„). For 
n = 2, it is known (probably at least since Sylvester [28]) that /(ai, 02) = 
aia2- For n > 2, all attempts for explicit formulas have proved elusive. 



2 



Here wc focus on the study of p[^{t), and show that it has an exphcit 
representation as a quasipolynomial. Through the discussion of we 
gain new insights into Probenius's problem. 

Another motivation to study p'j^{t) is the following trivial reduction 
formula to lower dimensions: 

P{ai,-..,an} W = J2 P'{a^,...,an-l}(t " "^«") " (3) 

m>0 

Here we use the convention that = if t < 0. This identity can be 

easily verified by viewing as 

= # |(mi,.. .,m„) e N" : ^m^afe = t-m„a„| . 

Hence, precise knowledge of the values of t for which = in lower 

dimensions sheds additional light on the Probenius number in higher di- 
mensions. 

The number appears in the lattice point count of 7^. It is for this 

reason that we decided to focus on this particular rational polytope. We 
present two methods (Sections 2 and 3) for computing the terms appearing 
in L{V°,t) and L{V,t). Both methods are refinements of concepts that 
were earlier introduced by the authors [2, 9]. In contrast to the mostly 
algberaic-geomctric and topological ways of computing L(V°, t) and L{V, t) 
[1, 6, 7, 14, 17, 18], our methods are analytic. In passing, we recover the 
Ehrhart-Macdonald reciprocity law relating L{V°,t) and L{V,t) [11, 20]. 
Within the computations a Dedekind-like finite Fourier sum appears, which 
shares some properties with its classical siblings, discussed in Section 4. In 
particular, we prove two reciprocity laws for these sums: a rederivation of 
the reciprocity law for Zagier's higher-dimensional Dedekind sums [30], and 
a new reciprocity law that generalizes a theorem of Gessel [13]. Finally, in 
Section 5 we give bounds on these generalized Dedekind sums and apply 
our results to give new bounds for the Frobenius number. The literature 
on such bounds is vast — see, for example, [4, 8, 12, 16, 25, 26, 27, 29]. 



2. THE RESIDUE METHOD 



In [2] , the first author used the residue theorem to count lattice points 
in a lattice polytope, that is, a polytope with integer vertices. Here we 
extend these methods to the case of rational vertices. 

We are interested in the number of lattice points in the tetrahedron V 
defined by (1) and integral dilates of it. We can interpret 

L{V,t) = # |(mi, . . . ,m„) G Z" : mfe > 0, ^rrifcafe < t| , 
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as the Taylor coefficient of z* of the function 

(1 + z"' + Z^"^ +...)••• (1 + + 0^"" + + z + z^ + ...) 

1 1 1 



1 - z^i I- z°-^l- z 
Equivalently 

Lip, t) = Res (- ^ — ^— — -, z = oV (4) 

If this expression counts the number of lattice points in tV, then the re- 
maining task is to compute the other residues of 

^-t-i 

F-t{z) 



(1 - z"^) • • • (1 - {I- z) 

and use the residue theorem for the sphere CU{oo}. F_t has poles at and 
all a*^, . . . , a}^ roots of unity. It is particularly easy to get precise formulas 
if the poles at the nontrivial roots of unity arc simple. For this reason, 
assume in the following that oi , . . . , a„ are pairwise relatively prime. Then 
the residues for the af^, . . . ,a^^ roots of unity are not hard to compute: 
Let A"i = 1 ^ A, then 

Res{F-t{z),z = A) = 

A-*-i „ / 1 



A-*-i z-\ 
lim 



(1 - A«2) ... (1 - A"") (1 - A) z-.\ 1 - 0«i 

A^^ 

ai(l-A°2)---(l- A»»)(l-A) ■ 

If we add up all the nontrivial a**^ roots of unity, we obtain 



Res{F_t{z),z = X) = 



«i A"±r^A (1-A«2)---(1-A«n)(i-A) 

c-kt 



«1 ^ (l-^fc«2)...(l_^fca„)(l_^fc) ' 

where ^ is a primitive a}^ root of unity. This motivates the following 
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Definition 1. Let ci, . . . ,c„ e Z be relatively prime to c e Z, and 
t G Ij. Define the Fourier-Dedekind sum as 



at{ci,...,Cn;c) = - y2 



A' 



c (A=i 



l)...(Ac.-i) 



Some properties of at are discussed in Section 4. With this notation, 
we can now write 

^ Res (F_((2), ^ = A) = (-1)"+V_t (a2, . . . , a„, 1; ai) . 

A<»i=l7iA 

We get similar residues for the d!^, • • • , a^'' roots of unity. Finally, note that 
Res(F_t, 2; = 00) = 0, so that the residue theorem allows us to rewrite (4): 

Theorem 1. Let V he given by (1), with ai, . . . ,a„ pairwise relatively 
prime. Then 

n 

L{V, t) = R-t (ai, . . . , a„) + (-1)" ^ cr_t(ai, . . . , fij, . . . , a„, 1; aj) 

where R-t (ai, . . . , fln) = —Res {F-t{z), z = 1), and aj means we omit the 
term aj. 

Remarks. 1. R-t can be easily calculated via 

Res iF-t{z), z = l)= Res {e''F-t{e^), z = 0) 



Res 



(1 - e«i^) ■••(!- e""^) (1 - 



,z = 



To facilitate the computation in higher dimensions, one can use mathe- 
matics software such as Maple or Mathematica. It is easy to see that 
R-t (ai, . . . , a„) is a polynomial in t whose coefficients are rational expres- 
sions in ai, . . . , a„. The first values for R-t are 



„ / ^ t 1 1 



t' 



t (\ 1 



a'^ b'^ ab 



-I- 



1 

12 



a _^ 6 _^ 1 \ 

b a ab J 

t' 



R-t{a,b,c) = —— + 



6abc 



t 



3 3 3 
12 Ka^b^c 



?( 


' 1 
ab 


1 

— -h 

ac 


1 

6^ + 


abc j 


3 


3 


3 


a 


b 


ab 


+ — 

ac 


^bc 




f — - 

ac 



1 1 

1 + - + - 

a b 



c _^ 1 

ab abc 



1 

24 



„333aa66ccl 1 1 a b c^ 
S+- + - + - + - + - + - + - + - + - + — + — + — + — + — + — 
a b c b c a c a b ab ac be be ac ab , 
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2. If ai, . . . , a„ arc not pairwise relatively prime, we can get similar formu- 
las for L(V,t). In this case we do not have only simple poles, so that the 
computation of the residues gets slightly more complicated. 

For the computation of L{V°,t) (the number of lattice points in the 
interior of our tetrahedron tV), we similarly write 

L{V°,t) = #|(mi,...,m„) e Z" : mfc > 0,^mfeafc < t| . 

So now we can interpret L{'P°,t) as the Taylor coefficient of of the 
function 

(^ai ^ ^2ax (^a„ ^ ^2a„ + _ _ /) + ^2 ^ _ _ ^ 

_ Z"' Z"^ Z 

~ 1 - 001 "' I- z"-n I- z ' 

or equivalently as 

Res z~*~^,z = 

= Res \ — ^ — ^ — - — z^^^.z = oo 

To be able to use the residue theorem, this time we have to consider the 
function 

The residues at the finite poles of Ft can be computed as before, with 
t replaced by —t, and the proof of the following theorem is completely 
analogous to Theorem 1: 

Theorem 2. Let V he given by (1), with ai, . . . , a„ pairwise relatively 
prime. Then 

n 

L{V°,t) = (-l)"i?t {au . . . , a„) + ^ crt(ai, . . . , a^, . . . , a„, 1; aj) 

As an immediate consequence we get the remarkable 
Corollary 1 (Ehrhart-Macdonald Reciprocity Law). 

L{V°,-t) = {-lYL{V,t) . 

This result was conjectured for convex rational polytopes by Ehrhart 
[11], and first proved by Macdonald [20]. 

Of particular interest is the number of lattice points on the boundary 
of tV. Besides computing L('P°,t) and L{V,t) and taking differences, we 
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can also adjust our method to this situation, cspcciaUy if wc arc interested 
in only parts of the boundary. As an example, we will compute as 
defined in the introduction (2), which appears in the context of the Frobe- 
nius problem. Again, for reasons of simplicity we assume in the following 
that ai , . . . a„ are pairwise coprime positive integers. 
This time we interpret 

PaW = #|(mi,...,m„) gN" : ^mfcttfe = i| 

as the Taylor coefficient of z* of the function 

{z"' + z''"' +...)■■■ (z""- + z"""- + .. .) 
_ 2;"! z"" 

~ 1 _ ^ai ■ ■ ■ 1 - ^"n ■ 

That is, 



p'^it) = Res 



-t-i 



= Res ( — 3_ 1 — z = oo 



Thus, we have to find the other residues of 



since 

p'^it) = - Res {Gt{z), z = oo). (5) 

Gt has its other poles at all af^, . . . , roots of unity. Again, note that Gt 
has simple poles at all the nontrivial roots of unity. Let A be a nontrivial 
af^ root of unity, then 

\t-i / 2 
Res (GAz), z = X) = — — — Res ( -, z = A 

^ ' (A''^ - 1) . . . (A"" - 1) V-^^'-l 

A* 



ai (A«2 -!)••• (A«" - 1) ■ 
Adding up all the nontrivial a^^ roots of unity, we obtain 

E Res(G.(.),. = A) = l (A^.-i).'.(A^„-i) 
= at (a2, . . . ,a„;ai) . 
Together with the similar residues at the other roots of unity, (5) gives us 
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Theorem 3. 

n 

PAit) ^ R't (ai,---,a„) + ^(7t(ai,...,aj, ...,a„;aj) , 

where R[ (ai, . . . , a„) = Res {Gt{z), z = 1). 

R' is as easily computed as before, the first values are 

R't{a,b,c) = 



b 

i / 1 1 1 



2abc 2 \ab ae be 
1/333 a b e 
12 V a 6 c be ae ab 



R[{a,b,c, d) 
+ 



f / 1 

Gabed 4 \ abe 
t/333333a b e d 



1 


+ ^ + M 

aed bed J 




abd 




3 


3 ^ a ^ 6 
ed bed aed 


^ e ^ 




abd 


b 

- — + 
ae 


bee 
ed ab ad 


e d 
bd ab 



12 \ab ae ad be bd ed bed aed abd abe 

aabbbeeeddd 

\ i.- + T3 + -i + —i + — + -i + -[: + —i + T3 + ^ + — + T- 

24 \bc bd ed ad ae ed ab ad bd ab ae be 
1111 

abed 

A general formula for R'^ (ai, . . . , a„) was recently discovered in [3]. 

For generalizations, note that we can apply our method to any tetrahe- 
dron given in the form (1), with the Ofe's replaced by any rational numbers. 
Moreover, any convex rational poly tope (that is, a convex polytope whose 
vertices have rational coordinates) can be described by a finite number of 
inequalities over the rationals. In other words, a convex lattice polytope 
V is an intersection of finitely many half-spaces. This description of the 
polytope leads to an integral in several complex variables, as discussed in 
[2, Theorem 8] for lattice polytopes. 



3. THE FOURIER METHOD 



In this section we outline a Fourier-analytic method that achieves the 
same results. Although the theory is a little harder, the method is of in- 
dependent interest. It draws connections to Brion's theorem on generating 
functions [5] and to the basic results of [9]. 

To be concrete, we illustrate the general case with the 2-dimcnsional ra- 
tional triangle P whose vertices are Wq = (0,0), t;i = (|,0), and = (O, |). 
As before, the number of lattice points in the 1-dimensional hypotenuse of 
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this right triangle is 

P{a,6}(*) = #{{m,n) gN^ :am + bn = t}. 

We denote the tangent cone to V at the vertex Vi by Ki. We recall that the 
exponential sum attached to the cone K (with vertex v) is by definition 

aK{s) = Yl e-'"<''™> , (6) 

where s is any complex vector that makes the infinite sum (6) converge. 
An equivalent formulation of (6) which appears more combinatorial is 

where = cc™^ • • • a;™" and xj = e~^'^*^' . 

In general dimension, let the vertices of the rational polytope V be 
vi, . . . ,vi. Let the corresponding tangent cone at vj be Kj . Finally, let the 
finite exponential sum over V be 

ap{s) = e-2-<^.-) . (8) 

Then there is the basic result that each exponential sum (7) is a rational 
function of x, and the following theorem relates these rational functions 
[5]: 

Theorem 4 (Brion). For a generic value of s € C", 

ar{s)=Y'^KM- (9) 

This result allows us to transfer the enumeration of lattice points in V 

to the enumeration of lattice points in the tangent cones Ki at the vertices 
of V, an easier task. In the theorem above, 'generic value of s' means any 
s e C" for which these rational functions do not blow up to infinity. 

To apply these results to our given rational triangle V, we first employ 
the methods of [9] to get an explicit formula for the exponential sum for 
each tangent cone of V. Then, by Brion's theorem on tangent cones, the 
sum of the three exponential sums attached to the tangent cones equals 
the exponential sum over V. Canceling the singularities arising from each 
tangent cone, and letting s ^ 1, we get the explicit formula of the previous 
section for the number of lattice points in the rational triangle V. 

In our case, Ki is generated by the two rational vectors —vi and V2—V1. 
We form the matrix 
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whose columns arc the vectors that generate the cone Ki . Once we compute 
o'Kiis), cfkq,{s) will follow by symmetry. The easiest exponential sum to 
compute is 



27r(misi+m2S2) 



m.,,>a 



(1 - e-^'^^i) (1 - e-^"*') " 

To compute (JKi{s) [i ^ 0), we first translate the cone Ki by the vector 
—Vi so that its new vertex is the origin. We therefore let K = Ki — Vi, 
and the following elementary lemma illustrates how a translation affects 
the Fourier transform. Let 



Xk{x) 



1 ifxeK, 
iix^K 



denote the characteristic function of K. 
Lemma 1. Let 

forxG M",s e C". Then 



Proof. 



Mo = I XK+v{x) e-^^^^'™) e2"<«'-> dx 
= f e2-'«+-'-> XK+v{x) dx 
= I e2-<«+-.^-"> XK{y) dy 



This lemma also shows why it is useful to study the Fourier transform 
of K at complex values of the variable; that is, at ^ + is. We study F{x) 
because (6) can be rewritten as 



10 



All of the lemmas of [9] remain true in this rational polytope context. The 
idea is to apply Poisson summation to X^^g^" ("^) ^^'^ write formally 

5^ F„(m) = ^ F„(m) 

The right-hand side diverges, though, and some smoothing completes the 
picture. Because the steps are identical to those in [9] , we omit the ensuing 
details. Let = . We get 



coth^ f si,i + ^ ) + 1 



where 



Si,i = (s, generator 1 of K^) = ( (si, S2) , ( --, 



a a 



and 

si,2 = (s, generator 2 of ATi) = ((si,S2) ,(--,,, , , , 

\ \ a J / a 

By (9), we have 

#{Z''ntV} = V 1 = limiaKois) + ctkAs) + (^kAs)) ■ 

mGZ^ntP 

Using the explicit description of OKi{s) in terms of cotangent functions, 
we can cancel their singularities at s = and simply add the holomorphic 
contributions to OKi (s) at s = 0. The left-hand side of (9) is holomorphic 
in s, so that we are guaranteed that the singularities on the right-hand side 
cancel each other. 

The only term in the finite sum (10) that contributes a singularity at 
s = is the r = term. We expand the three exponential sums (iKiis) 
into their Laurent expansions about s = 0. Here we only require the first 
3 terms of their Laurent expansions. In dimension n we would require the 
first n-\-\ terms; otherwise every step is the same in general dimension n. 

We make use of the Laurent series 

near s = 0, as well as the Laurent series for cotTrs near s = 0. After 
expanding each cotangent in (10) for gkq {s), (TKi (s) and (s) and letting 
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s —>■ 0, we obtain Theorem 1 above as 



1 / 1 1\ 1 /a b 1 
+ + - + - - + - + — 

4\ a b J 12 \b a ab 

1 ^ 1 f 



Note that, as before, the perfodic portfon of L{V, t) is entirely contained 
in the "constant" t term. By Ehrhart's reciprocity law (Corollary 1, [11]), 
there is a similar expression for L{V,t), and taking 





't' 




't' 






a 




b 





1 



L(V,t)-L{V,t) 
gives us P{o,6}(i)- The same analysis gives us Theorem 1 in M". 

4. THE FOURIER-DEDEKIND SUM 

In the derivation of the various lattice count formulas, we naturally 
arrived at the Fourier-Dedekind sum 

a,(c„...,c„;c) = l (AC - 1) (A- - 1) " 

This expression is a generalization of the classical Dedekind sum s(/i, fc) 
[23] and its various generalizations [10, 13, 21, 22, 30]. In fact, an easy 
calculation shows 



ao(a,l;c) = l ^ I = 



c-1 



(A- - 1) (A - 1) 

111 ^ Trfca tt/c 1 1 , , 

- - — > cot cot — = - - 5(a, c) . 

4 4c 4c ^ c c 4 4c 

fe=i 

In general, note that ctj (ci, . . . , c„; c) is a rational number: It is an element 
of the cyclotomic field of c*'^ roots of unity, and invariant under all Galois 
transformations of this field. 
Some obvious properties are 

CTt (ci, . . . , c„; c) = at (c^(i) , . . . , c^(„) ; c) for any tt G 5„ 

ut (ci, . . . , c„; c) = (T(t mod c) (ci mod c, . . . , c„ mod c; c) (11) 

CTt (ci, . . . , c„; c) = CTfet (6ci, . . . , 6c„; c) for any 6 e Z with (6, c) = 1 
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We can get more familiar-looking formulas for at in certain dimensions. 
For example, counting points in dimension 1, we find that 



L{V, t) = # {m G Z : m > 0,mc < t} = 
so that Theorem 1 implies 



i7-t(l;c) 



1 



E 



A-* t 


t 


1 


(A-1) - c 


_c_ 


2 



+ 



2c 



2c 



• (12) 



Here, {{x)) — x — \ x\ — 1/2 is a sawtooth function (differing slightly from 
the one appearing in the classical Dedekind sums). This restates the well- 
known finite Fourier expansion of the sawtooth function (see, e.g., [23]). 
As another example, we reformulate 



{a,h;c) = - ^ 



by means of finite Fourier series. Consider 

c-l 

I — , A^" I 

<Jt{a;c) 



1 ^ ^ka-'^t 



1 

2c 



(13) 



where ^ is a primitive c**^ root of unity and aa~^ = 1 mod c; here, the last 
equality follows from (12). We use the well-known convolution theorem for 
finite Fourier series: 

Theorem 5. Let f{t) = ^Ef=o'«fe^'* (^^d g{t) = jjEk^ohe*, 
where ^ is a primitive AT*'* root of unity. Then 



N-1 



Hence by (13), 



fe=o 



c-l 



m=0 



at{a, b;c)= ^ crt-m{a; c)am{b; c) 



E 

c-l 

E 

m=0 



—a ^{t — m) 



a ^{m — t) 



1 

^2c 



1 

4c 



1 ■ 
2^ 
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Here, aa ^ = bb ^ = 1 mod c. The last equality follows from 

tm-i 

m=0 

Furthermore, by the periodicity of ((a;)), 

m=Q ^ ^ ' ' 

The expression on the right is, up to a trivial term, a special case of a 
Dedekind-Rademacher sum [10, 19, 21, 22]. It is a curious fact that the 
function (Tt(a, b; c) is the nontrivial part of a multiplier system of a weight- 
modular form [24, p. 121]. 

We conlude this section by proving two reciprocity laws for Fourier- 
Dedekind sums. The first one is equivalent to Zagier's reciprocity law for 
his higher dimensional Dedekind sums [30]. They are essentially Fourier- 
Dedekind sums with t = 0, that is, trivial numerators. 

Theorem 6. For pairwise relatively prime integers ai, . . . , o„, 

n 

cro(ai, . . . , Oj, . . . , a„; aj) = 1 — (ai, . . . , a„) , 

i=i 

where R[ is the rational function given in Theorem 3. 

Proof. It is well known [11] that the constant term of a lattice polytope 
(that is, a polytope with integral vertices) equals the Euler characteristic 
of the polytope. Consider the polytope 

|(a;i,...,x„) e M>o : ^Xfeafe = l| , 

whose dilates correspond to the quantor of Theorem 3. If we dilate 

this polytope only by multiples of ai • • • a„, say t = ai • • • a„w, we obtain 
the dilates of a lattice polytope. Theorem 3 simplifies for these t to 

n 

P'a{(^1 ■ ■ ■ dnW) = R'a^.-anw (oi , • • • , «n) + X] ^o(ai, ■ ■ ■ , CLj , . . . , an, ttj) , 

using the periodicity of at (H)- On the other hand, wc know that the 
constant term (in terms of w) is the Euler characteristic of the polytope 
and hence equals 1, which yields the identity 

n 

1 = i?o (ai, • • • ,an) + X^o-o(ai, ...,aj,.. . ,a„;aj) . 
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The second one is a new reciprocity law, which generalizes the following 

[13] 

Theorem 7 (Gessel) . Let m and n be relatively prime and suppose that 
< r < m + n. Then 



1 V A''+^ 1 ^ 

m_ ^ f A" - 1) (A - 1) n , ^ 



- - 1) (A - 1) n ^J^^^ (A™ - 1) (A - 1) 

- + i-l 



h 


11 




n 


m 


mn ) 


(1 


1 

+ - 




V m 


n 


mn 7 



It is not hard to see that Gessel's theorem follows as the two-dimensional 
case of 

Theorem 8. Let ai,...,a„ he pairwise relatively prime integers and 
<t <ai-\ h a„. Then 

n 

^ crt(ai, . . . , dj", . . . , a„; Uj) = —R[ (ai, . . . , a„) , 

where R[ is the rational function given in Theorem 3. 

Proof. By definition, = if < f < ai H hfln- Hence Theorem 

3 yields an identiy for these values of t: 

n 

= Rf (ai, ...,a„) + ^crt(ai, ...,aj,.. .,an;aj) . 



It is worth noticing that both Theorems 6 and 7 imply the reciprocity 
law for the classical Dcdckind sum s(a, 6). It should be finally mentioned 
that in special cases there arc other reciprocity laws, for example, for the 
sum appearing on the right-hand side in (14) [10, 22]. Wc note that, as a 
consequence, we can compute (T((a, b; c) in polynomial time. 



5. THE FROBENIUS PROBLEM 



In this last section wc apply Theorem 3 (the explicit formula for p^(t)) 
to Frobcnius's original problem. As an example, we will discuss the 3- 
dimensional case. Note that a bound for dimension 3 yields a bound for 
the general case: It can be easily verified that 

/(ai, . . . ,a„) < /(ai,(X2,a3) + a4 H l-a„ (15) 
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Furthermore, in dimension 3 it suffices to assume that ai, a2, as are pairwise 
coprime, due to Johnson's formula [15]: If <? = (ai, 02), then 



/(ai, 02,03) =3- / ( — ,— ,03 
9 9 



(16) 



Now assume a, b, c pairwise relatively prime, and recall (14): 



at{a,b;c) = ^ 



m=0 



-a-'^{bm + t) 



1 

4c 



where aa = 1 mod c. We will use the Cauchy-Schwartz inequality 



(k) 



k=l 



< 



(17) 



fe=i 



Here Ofe G M, and tt G iS^ is a permutation. Since (a ^6, c) = 1, we can use 
(17) to obtain 



TO=0 



m=0 



m 1 

7^2 



1 

4^ 



1 (2c- l)(c- l)c lc(c-l)_c_l 
6 c 2 4 ~ 4c 

c 1 

12 ~ 12^; ■ 



This also restates Rademacher's bound on the classical Dedekind sums [23]. 
Using this in the formula for dimension 3 (remark after Theorem 3), we 
get 



P'{a,b,c}(t) > 



2a6c 2 \ ah ac be 



1/333 a b c 
12 \ o b c be ac ab 



-J^ia + b + c) 



1/111 

12 \a^b^c 
1 fa 



t / 1 ^ 1 ^ 1 
2a6c 2 \ab ac be J 12 \bc ' ac ab 

1 / . X 1 A 1 1 
12 6 \ a c 



The larger zero of the right-hand side is an upper bound for the solution 
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of the Frobenius problem: 



/(a, b, c) < abc 



1 / 1 



+ 



1 



1 



2 V o6 be ac 



1/1 1 1 
4 \ab be ae 



1 



12 







1 










-b' 


-I) 


)] 



1/; 



< - (a + 6 + c) + afec-i 



'i f 1 1 1 
4 \ab be ac 



1/1 1 1 
6 \ab be ac 



1 / 1 



{a + b + e)+abeJ- {- + - + - 



1 



< - (a + 6 + c) + a6c4 



2 \a& 6c ac 
4 V ao oc ac 



1 f — + - + — ^ + - 

2 \ a6 6c ac y 3 



For the last inequality, we used the fact that ^ + ^ + ^ < ^'^To'^TE ~ 5' 
This proves, using (15) and (16), 

Theorem 9. Let a\ < a2 < ■ ■ ■ < an be relatively prime. Then 

f{ai, . . . , a„) < ^ ^VaIa^a7(a7Ta2"-i-a3) + oi + 02 + asj +04 H |-a„ . 

Remarks. 1. Sometimes the Frobenius problc;m is stated in a slightly 
different form: Given relatively prime positive integers ai, . . . , a„, find the 
largest value of t such that X)fe=i "^fc'^fe = * solution in nonnega- 

tive integers mi, ... , m„. This number is denoted by g{ai, . . . , a„). It is, 
however, easy to see that 

g{ai, . . . , a„) = f{ai, . . . , a„) - ai a„ . 

So we can restate Theorem 9 in a more compact form as 

g{ai, . . . , a„) < i (^\/ 010203 (ai + 02 + as) - oi - 02 - 03^ . 

2. Bounds on the Frobenius number in the literature include results by 
Erdos and Graham [12] 



Selmer [27] 
and Vitek [29] 



g(oi, . . . ,o„) < 2a„ 
/(oi, . . . ,a„) < 2a„_i 



Oi 



On 

n 



g{ai, ...,a„)< 



-(o2-l)(a„-2) 



1 . 
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Theorem 9 is certainly of the same order. What might be more interesting, 
however, is the fact that the bound in Theorem 9 is of a different nature 
than the bounds stated above: namely, it involves three variables, and is 
thus — especially in terms of estimating g{ai, 02, as) — more symmetric. 
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